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S-H ■ Abstract. We present evidence for a close analogy between the nonlinear behavior of a pulsed 

i-rj ' microwave-driven Josephson junction at low temperature and the experimentally observed be- 

havior of Josephson systems operated below the quantum transition temperature under similar 
conditions. We specifically address observations of Ramsey-type fringe oscillations, which can be 
understood in classical nonlinear dynamics as results of slow transient oscillations in a pulsed 
microwave environment. Simulations are conducted to mimic experimental measurements by 

*vj ' recording the statistics of microwave-induced escape events from the anharmonic potential well of 

a zero- voltage state. Observations consistent with experimentally found Ramsey-type oscillations 

— , ■ are found in the classical model. 

^h ; 
^- 

O ' 1. Introduction 

Macroscopic quantum behavior of Josephson systems has been the focus of 
intense research over the past years due to the potential applications of Josephson 
technology in quantum information processing [1]. Many experimental observa- 
tions have been conducted in order to characterize and understand the quantum 
properties of various Josephson device and networks [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. 
Due to the nature of the Josephson effect it is not easy to make direct observations 
of the expected quantum states in the anharmonic potential of a zero- voltage state 
Josephson system, so most measurements are conducted statistically by recording 
the distributions of transitions from one system state to another. These transi- 
tions are commonly induced through the controlled application of microwaves and 
microwave pulses of frequencies that are commensurate with the expected energy 
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separations between possible quantum levels of the system. 

Some of the key microwave induced features that have been used to illustrate 
the quantum properties of Josephson systems are multi-peaked switching distribu- 
tions as a function of bias-point in continuously applied microwaves [13, 14], obser- 
vations of Rabi-oscillations [15] in the switching statistics of systems perturbed by 
pulsed microwaves [3, 7, 12], and the associated Ramsey- fringes [16] when several 
coordinated and parameterized microwave pulsed are applied sequentially [10, 11]. 
These observations, along with many others, have produced a large body of in- 
sight to the dynamical response of Josephson systems to microwave perturbations 
applied at very low temperature, and the acquired data has been used to inter- 
pret the quantum nature of the superconducting networks under investigation. 
However, many of the observations may also be attributed to classical nonlinear 
dynamics and statistics. For example, it has been demonstrated that the multi- 
peak switching distributions as a function of varying system bias for a constantly 
applied microwave field can be understood as classical high-Q resonances in the 
nonlinear system [14, 17]. Also Rabi-type oscillations in the switching statistics of 
Josephson systems perturbed by microwave pulses have been shown to be consis- 
tent with theoretical and numerical expectations based on a classical model of a 
Josephson junction, which exhibit the transients in response to the pulsed pertur- 
bations [18, 19]. It is the purpose of this paper to demonstrate that Ramsey-type 
fringes can also be observed in the classical model when the recipe of the exper- 
imental measurements is followed. We here, for simplicity, illustrate the results 
with the simplest possible system, namely a single classical Josephson junction, 
perturbed by a bias-current, microwave fields, probe fields, and thermal noise. 

2. The Model 

The normalized classical equation for a perturbed Josephson junction within 
the RCSJ model (see Figure la) is [20], 

ip + aip + sin (p = r] + e s (t)sm(u! s t + S ) + e p (t) +n(t) , (1) 

where ip is the difference between the phases of the quantum mechanical wave 
functions defining the junction, r\ represents the dc bias current, and e s (t), ui s , and 
9 S represent microwave amplitude, frequency, and phase, respectively (see Figure 
lb). All currents are normalized to the critical current I c , and time is measured in 
units of the inverse plasma frequency, wo -1 , where ujq 2 = 2eI c /fiC = 2itI c /&oC, 
C being the capacitance of the junction, and <&o = h/2e is the flux quantum. 
Tunneling of quasiparticles is represented by the dissipative term, where a = 
htuo/2eRI c is given by the shunt resistance R, and the accompanying thermal 
fluctuations are defined by the dissipation-fluctuation relationship [21] 

(n(t)) = (2) 
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Fig. 1: (a) Sketch of the RCSJ model where I/I c = rj + e s (t) sin(uj s t + 9) + e p (t) 
and N = n(t)I c . (b) Sketch of the current supply to the junction. 
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T being the thermodynamic temperature, Hj is the characteristic Josephson en- 
ergy Hj = I c h/2e, and is thereby defined as the normalized temperature. A 
current pulse for probing the state of the system is represented by e P (t). 



3. Simulation Details 

Following the procedure of reported experiments, we record the switching from 
the zero- voltage state ((0) = 0) as a result of applying the external current 
sketched in Figure lb. We first equilibrate the system at a chosen temperature for 
a given value of rj. For a randomly chosen, but temporally constant, phase 9 S , a 
microwave pulse with frequency uj s , which is associated with the natural resonance 
of the junction at bias point rj, is applied for a duration of At a- This duration is 
chosen to be At a = Tr/4, where tr = 2-k/Qr is the Rabi-type oscillation period 
[18, 19] for the microwave amplitude e s . This pulse is denoted a 7r/2-pulse due 
to its parameterization based on Rabi-type oscillations. Another 7r/2-pulse, in 
phase 9 S with the first pulse, is applied at a later time with identical amplitude, 
frequency, and duration (A£g = At a), whereafter a short pulse e p (i) is applied 
to probe whether the system is in a high or low energy state after the sequential 
microwave pulse application. The probe pulse is parameterized such that a state 
of relatively high system energy at £3 will result in a likely escape event when e p (t) 
is activated, while a state of relatively low system energy will result in a low prob- 
ability of escape when e p (t) is activated. This procedure is repeated many times, 
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Fig. 2: Josephson response to the sequential application of two ir/2 microwave 
pulses followed by a probe field. (a,c) Josephson phase. (b,d) System energy as 
defined by Equations (4) and (5). Parameters are a = 1CP 4 , r] = 0.904706, e s = 
2.17- 10" 3 , e p = 8.2 -10~ 2 , and = 0. (a,b) Non-switching for At d = t 2 -ti = 400. 
(c,d) Switching for Atd = ^2 — t\ = 70. 

each time with a new realization of 9 S , in order to generate information about 
the probability of switching for a given set of parameters, such as the time delay 
between the two 7r/2-pulses. Figure 2 illustrates typical examples of the events de- 
scribed above for = 0. Specific parameter values are: a = 10 -4 , e s = 2.17- 10~ 3 , 
e p = 8.2 • 10 -2 , lo s = uji, and rj = 0.904706, which results in a linear resonance 
frequency of u\ = \/l — r/ 2 = 0.652714. 7r/2-pulses are applied in the intervals 
t E [0;*i = At A ] and t € [t 2 ;t 3 = t 2 + At B ] with At A = At B = 90.75 time units. 
This value is obtained from the corresponding Rabi-type oscillations reported in 
[19]. After the second 7r/2-pulse, a short delay of At = 20 is allowed before the 
probe pulse £ p (t) is initiated at £4. The probe has linear rise and fall times of 160 
normalized time units with a constant value interval of 275 time units. We show 
the system response to the microwave perturbations in terms of both phase (p and 
the normalized energy H — Hq, which is defined by 



H = -ip + 1 — cos ip — r/cp 



Ho 



•q 2 — rj sin 7] . 



(4) 
(5) 
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Fig. 3: Switching probability as a function of 7r/2-pulse separation for = 2- 10~ 4 . 
Each point represents statistics of 2,500 events. Parameters are as in Figure 2, 
unless otherwise listed, (a) a = 10~ 4 . Visible Ramsey-type fringe oscillation 
for lj s = uji = \/l — f] 2 = 0.652714. Inset shows vanishing Ramsey-type fringe 
frequency near uj s = 0.990a;/. (b) same as (a), except for a = 10 -3 . Inset shows 
vanishing Ramsey-type fringe frequency near to s = 0.997a;;. 



Figures 2a and 2b show the response for Atd = 400. It is evident that the first 
7r/2-pulse elevates the system energy by leaving the system in a phase-locked state 
at t = £2 with energy E(t\). At the onset of the second 7r/2-pulse, a significant 
phase-slip between the oscillation of <p and the microwave has developed during 
the interval [£1; £2], and the second microwave pulse therefore decreases the energy 
of the Josepshon system. The system is left at t = £3 with energy E{t^) < E(ti), 
which, when the probe field is applied at £4, results in only temporary energy 
increase while the probe is applied; this is not enough to make the system switch 
into a non-zero voltage state, which can be seen from the fact that E{t) — > 
for large £). For comparison we show Figures 2c and d for A£^ = 70. We here 
observe the same initial behavior, but the application of the second 7r/2-pulse leaves 
the system in a relatively high energy state E{t^) > E(t\), since the phase-slip 
between the junction and the microwave field is minor in this case. Consequently, 
the junction phase cp switches out of the bound state when the probe pulse is 
applied at £4. This is recognized by the diverging energy for £ > £4. 
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Fig. 4: Ramsey-type fringe frequency Qp as a function of applied microwave fre- 
quency for two different dissipation parameters: (a) a = 10~ 4 , and (b) a = 10~ 3 . 
Other parameters are as in Figure (3). Arrows indicate the driving frequency for 
which lu s = {jj[. Dashed lines indicate a slope of ±1. 



4. Ramsey- Type Fringes 

Ramsey-type fringes in the switching probability as a function of 7r/2-pulse 
separation can be directly observed in Figure 3, where we have shown the variation 
of the switching probability for = 2- 10 -4 , calculated as averages of 2,500 events 
for different values of a and u s . Figure 3a clearly shows a distinct frequency, which 
we will name the Ramsey- type fringe frequency dp, for uj s = uj\ 



<n 



T 



0.652714 and a = 10 -4 . The inset shows that the fringe frequency depends on the 
applied microwave frequency such that to s = 0.990a;; results in £lp ~ 0. Similar 
behavior is observed for larger dissipation parameter a = 10~ 3 shown in Figure 
3b. Clearly, the fringe-frequency is different for uj s = u>i, compared to Figure 3a, 
and in this case to s = 0.997a;; results in Qp ~ 0. 

A number of simulations of the Ramsey-type fringe frequency for different 
microwave frequencies have been conducted, and the results are summarized in 
Figure 4, where the fringe frequency 0,f is shown as a function of microwave 
frequency lo s for the two different dissipation parameters mentioned in Figure 3. 
The vertical arrows indicate the linear resonance frequency uj\ = \J\ — rj 2 . The 
simulation data are represented by markers with error bars. The frequencies are 
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read from figures like Figure 3. For low £lp, less than one wavelength of fringe 
oscillation may be visible and the oscillation may be somewhat distorted (see 
Figure 3b). In such cases we measured the time from the first trough (peak) to 
the next peak (trough) and used the half period to determine the frequency. In the 
case of saturation such as Figure 3a, we measured the period of the first regular 
waveform. These issues also account for some of the distortion in Figure 4. 

We observe the distinct "V-shape" signature of the Ramsey- type frequency 
Qf(u> s ) [10] with a slope very close to ±unity (lines with slope ±1 are shown in 
Figure 4 along with the simulation data), and tip ~ for microwave frequencies 
close to wj. It is noticeable that the characteristic microwave frequency for which 
Qf = increases with a, and we further notice that this characteristic frequency 
is smaller than u\. While we have found no experimental data for comparison, 
future experiments may want to address this observation. 

5. Conclusions 

Like Rabi-type oscillations and resonant multi-peak switching distributions, 
Ramsey-type fringes have been used to identify and analyze expected macroscopic 
quantum behavior of Josephson systems. In this paper we have addressed the com- 
parable classical system through direct simulations of the well-established nonlin- 
ear classical model equation, which is driven according to the recipe prescribed by 
the experimental reports of Ramsey- fringes. Our results show that even the sim- 
plest possible classical Josephson model, the single RCSJ model, clearly exhibits 
the features of Ramsey-type fringes when the model is used to simulate a low 
temperature, low dissipation system. A detailed theoretical analysis of classical 
Ramsey-type fringes [22] can be developed in close analogy with, and extention of, 
the analysis of phase-locking transients that have been found to be responsible for, 
e.g., Rabi-type oscillations in pulsed microwave-driven classical Josepshon systems 
[18, 19]. 
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